The Groverian Measure of Entanglement for Mixed States 
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The Groverian entanglement measure introduced earlier for pure quantum states [O. Biham, 
M.A. Nielsen and T. Osborne, Phys. Rev. A 65, 062312 (2002)] is generalized to the case of 
mixed states, in a way that maintains its operational interpretation. The Groverian measure of 
a mixed state of n qubits is obtained by a purification procedure into a pure state of 2n qubits, 
followed by an optimization process before the resulting state is fed into Grover's search algorithm. 
The Groverian measure, expressed in terms of the maximal success probability of the algorithm, 
provides an operational measure of entanglement of both pure and mixed quantum states of multiple 
qubits. These results may provide further insight into the role of entanglement in making quantum 
algorithms powerful. 

PACS numbers: 03.67.Lx, 89.70.+C 



I. INTRODUCTION 

The potential speedup offered by quantum computers 
is exemplified by Shor's factoring algorithm [l|, Grover's 
search algorithm @,0> and algorithms for quantum simu- 
lation . Although the origin of this speed-up is not fully 
understood, there are indications that quantum entangle- 
ment plays a crucial role in making quantum algorithms 
efficient In particular, it was shown that quantum 

algorithms that do not create entanglement can be sim- 
ulated efficiently on a classical computer 0. It is there- 
fore of interest to quantify the entanglement produced 
by quantum algorithms and examine its correlation with 
their efficiency. This requires to develop entanglement 
measures for the quantum states of multiple qubits that 
appear in quantum algorithms. These include pure states 
as well as mixed states, which would inevitably show up 
when decoherence effects are taken into account. 

The special case of bi-partite entanglement has been 
studied extensively in recent years and suitable entangle- 
ment measures were introduced. It was established that 
bi-partite entanglement can be considered as a resource 
for teleportation 8] . The entanglement of bi-partite pure 
states can be evaluated by the von Neumann entropy of 
the reduced density matrix, traced over one of the par- 
ties. For bi- par tite mixed states, several measures were 
proposed [2, flOl Hil l and for the special case of states of 
two qubits an exact formula for the entanglement of for- 
mation was obtained |121 1 1 3| . For mixed states of mul- 
tiple qubits, entanglement measures based on distance 
measures in Hilbert space were proposed EH El • 

Consider a mixed quantum state p of n qubits. The 
state is non-entangled, or separable, if its density matrix 
can be written in the form 



(i) 



where p„, k = 1, . . . , n is a density operator of a pure 
state of the fcth qubits, namely = and 
5^ P M = 1. In the special case that p is a pure state, all 



probabilities vanish except for Pi = 1, and the state can 
be expressed by 



= \il> 1 )®---®\il> n ). 



(2) 



Such states are called tensor-product states. In or- 
der to evaluate the entanglement of a quantum state, 
p, one needs a scalar function E(p) [or E(ip) for pure 
states] called an entanglement measure that satisfies 
EH E [H E3: (a) E(p) = if and only if p is 
a separable state; (b) Assuming that each qubit is held 
by a different party, it is not possible to increase E{p) 
by local operations and classical communication (LOCC) 
between the parties. Consider the special case of local 
unitary operators. Such operators cannot decrease E(p) 
because if they could then the inverse operators (which 
are also unitary) would increase it and thus contradict 
the second condition above. The conclusion is that local 
unitary operators cannot change E(p). 

The Groverian entanglement measure, G(tp), provides 
an operational measure of entanglement for pure states 
of multiple qubits EH • It is related to the success prob- 
ability of Grover's search algotirhm when the state \ip) 
is used as the initial state. A pre-precessing stage is al- 
lowed in which an arbitrary local unitary operator is ap- 
plied to each qubit. These operators are optimized in 
order to obtain the maximal success probability of the 
algorithm, P max (ip). The Groverian measure is given by 
G(tp) = y/l — P ma x(tp)- The Groverian measure was 
used in order to evaluate the entanglement in certain 
quantum states of high symmetry as well as in states 
that are gene rated during the operation of quantum al- 
gorithms |2Cj . For example, it was found that Grover's 
iterations generate highly entangled intermediate states, 
even in case that the initial and the final states are prod- 
uct states. 

In this paper we generalize the Groverian entanglement 
measure to the case of mixed states. The Groverian mea- 
sure, G(p), of a given mixed state p, of n qubits, is ob- 
tained by its purification into a pure state of 2n qubits. 
An optimization procedure based on Uhlmann's theorem 
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[2l| is then applied before the resulting pure state is fed 
into Grover's algorithm. G(p) is then expressed in terms 
of the maximal success probability P m ax(p), as described 
above for pure states. 

In Sec. [n] we briefly describe Grover's search algo- 
rithm, in a context suitable to this paper. In Sec. IIIII 
we review the Groverian entanglement measure for pure 
states. The generalization to mixed states is presented in 
Sec. IIVI and its operational interpretation is considered. 
The results are summarized and discussed in Sec. [3 



II. GROVER'S SEARCH ALGORITHM 



Using Eq. (J3J we obtain 

P 8 (i>) = \(r 1 \ij)\ 2 + 0(1/N), (6) 

namely, the success probability is determined by the over- 
lap between \tp) and |ry). 

III. THE GROVERIAN MEASURE FOR PURE 
STATES 



Consider a search space D containing N elements. We 
assume, for convenience, that N = 2", where n is an in- 
teger. This way, the elements of D can be represented by 
an n-qubit register \x) — \xi,X2, ■ ■ ■ ,x n ), with the com- 
putational basis states i = 0, . . . , N — 1. We assume 
that one element in the search space is marked, namely 
it is the solution of the search problem. The distinction 
between the marked and unmarked elements is expressed 
by a suitable function, / : D — ► {0, 1}, such that / = 1 
for the marked element, and / = for the rest. The 
search for the marked element now becomes a search for 
the element for which / = 1. To solve this problem on a 
classical computer one needs to evaluate / for each ele- 
ment, one by one, until a marked state is found. Thus, 
on average, N/2 evaluations of / are required and N 
in the worst case. For a quantum computer, on which 
/ is evaluated coherently, it was shown that a sequence 
of unitary operations, called Grover's algorithm and de- 
noted by {Jo, can locate the marked element using only 
0{\fN) coherent queries of / 0,0- 

Starting with the equal superposition state, 



\V) 



1 N - 1 



(3) 



and applying the operator Uq one obtains Ug\v) = 
\m) + 0(1/N), where \m) is the marked state. Thus, 
the success probability of the algorithm is almost unity 
0, 13- With this performance, Grover's algorithm was 
shown to be optimal |22l namely, it is as efficient as the- 
oretically possible [23. The adjoint equation takes the 
form 



(ry| = (m\U G + 0(l/N), 



(4) 



where the error is due to the discreteness of the Grover 
iterations 0. If an arbitrary pure state, is used 
as the initial state instead of the state \rj), the success 
probability is reduced |2^, |2(j . It is given by |53 



Consider Grover's search algorithm, in which an arbi- 
trary pure state \ip) is used as the initial state. Before 
applying the operator Uq, there is a pre-processing stage 
in which arbitrary local unitary operators, U\, U2, ■ ■ . , 
U n , are applied on the n qubits in the register (Fig. 
These operators are chosen such that the success proba- 
bility of the algorithm would be maximized. The maxi- 
mal success probability is thus given by 



PmaxW>)= max \(m\U G (U 1 ^---^U n )\^}\ 2 . (7) 
Ui , U2 , ■ ■ ■ , u n 



Using Eq. Q, this can be re- written as 



PmaxWO = max KvPl® 

Ul,U 2 ,...,Un 



u n m\ (8) 



or P maK (ip) — max|0) eT K^l?/')! 2 , where T is the space of 
all tensor product states of the form \cj)) = |0i)<8>- • -<8>|</>n)- 
The Groverian measure is given by [19j 



G(i/>) 



1 — max KqMiMI 2 . 



(9) 



P a W) = \(m\U G m | 2 +Q(1/N). 



(5) 



For the case of pure states, for which G(ip) is defined, it is 
closely related to an entanglement measure introduced in 
Refs. PHG3 and was shown to be an entanglement 
monotone. The latter measure is defined for both pure 
and mixed states. It can be interpreted as the distance 
between the given state and the nearest separable state 
and expressed in terms of the fidelity of the two states. 
Based on these results, it was shown [Hjj that G(ip) sat- 
isfies: (a) G(ip) > 0, with equality only when is a 
product state; (b) G(ip) cannot be increased using local 
operations and classical communication (LOCC). There- 
fore, G(tp) is an entanglement monotone for pure states. 
A related result was obtained in Ref. [2^|, where it was 
shown that the evolution of the quantum state during 
the iteration of Grover's algorithm corresponds to the 
shortest path in Hilbert space using a suitable metric. 
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FIG. 1: The quantum circuit that exemplifies the operational 
meaning of G(ip). An arbitrary pure state \ip) of n qubits is 
inserted as the input state. In the pre-processing stage, a 
local unitary operator is applied to each qubit before the re- 
sulting state is fed into Grover's algorithm. The local unitary 
operators Ui, i = 1, . . . , n are optimized in order to maximize 
the success probability of the search algorithm for the given 
initial state \ib). 



maximal success probability is given by 



-Pmax(p) = maxmaxmax|(m|[/GC^<A|^)| 



Using Eq. J3J we obtain 



max max max | (77 1 \ ip) \ A 



(10) 



(11) 



where S is the set of separable states of n qubits. The 
maximization is over all separable states a of n qubits, 
and for each of them, over all possible purifications \<p) 
of 2n qubits. This can be rewritten as 



-Pmax(p) = max max I (<^) I ^ 

\<t>) m 



(12) 



The first maximization is over all possible states \4>) of 
2n qubits which are purifications of separable states a of 
n qubits. The second maximization is over all states \tp) 
of 2n qubits which are purifications of p. According to 
Uhlmann's theorem, the fidelity of any two states p and 
a of ri qubits satisfies 



IV. THE GROVERIAN MEASURE FOR MIXED 
STATES 

A. Desired Properties and construction of G(p) 

The quantum circuit that demonstrates the evaluation 
of the Groverian measure G(p) for a mixed state p is 
shown in Fig. [5] The given state p is purified into a pure 
state of 2n qubits. In the pre-processing an operator 
Utf, is applied. This operator has the property that its 
adjoint operator satisfies C/j|?y) = \4>), where \(f>) is a pu- 
rification of a separable state of n qubits. The resulting 
state is then inserted into Grover's search algorithm. An 
optimization is performed over the entire family of op- 
erators U4, that satisfy the above condition, in order to 
maximize the success probablity of the search algorithm. 
Equivalently, the optimization can be performed over all 
the separable states, a, and over all possible purifications 
\cf>) of each one of them. The maximal success probabil- 
ity for the initial state p is denoted by -P m ax(p) and the 
Groverian measure is given by G{p) = yl — P m a,x(p)- 

The construction of the Groverian measure for mixed 
states is based on a purification process, where the pure 
state l^) of 2n qubits is a purification of the mixed state 
p of n qubits. Similarly, \4>) is a purification of a sepa- 
rable mixed state o. As in the case of pure states, we 
introduce a pre-processing stage before the state is in- 
serted into Grover's algorithm. During pre-processing an 
optimization is performed over a certain class of unitary 
operators U^. These operators satisfy \(j>) — U^\i]), where 
\4>) is a purification of a separable state a of n qubits. The 



F(p, a) = max max |{^>|^) | 2 , (13) 
\<t>) W 

where \4>) and \ip), of 2n qubits, are purifications of p 
and cr respectively j^J. A useful corollary (presented in 
Excercise 9.15 on page 411 of Ref. 3) enables us to 
remove the optimization on \tp), leading to 

F(p,a) = max | («^) | 2 , (14) 
\<t>) 

where is an arbitrary purification of p. Using this 
corollary we find that 

G(p) = ^l-maxF(p,a). (15) 

An entanglement measure for mixed states, p, should sat- 
isfy the conditions: (a) G(p) > 0, where the equality is 
obtained only if p is a separable state; (b) G(p) cannot 
be increased using LOCC. To express the second condi- 
tion we introduce a complete set of operators, defined by 
{^i}™ 1 where M, = M t {\) <g> Mi{2) <g> • • • <g> Mi{n) and 
YlT=i MiM^ = I. In this notation, the condition is that 
for every density matrix p 



G\Y^M iP M}\<G(p). (16) 

We proceed to prove that G(p) fulfills the two condi- 
tions described above: As shown in Refs. 0, 0, Il6| . 
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takes the form 



Purify 



U<t> 



u 



Measure 



FIG. 2: The quantum circuit for the Groverian measure G(p) 
of a mixed state of n qubits. The given state p is purified into 
a state |^>) of 2n qubits. In the pre-processing stage, a unitary 
operator U$ is applied on \ip), before the resulting state is fed 
into Grover's search algorithm in the space of 2n qubits. The 
operator belongs to a class of operators that their adjoint 
operators satisfy U\\r}) = \<j>) where \<j>) is a purification of a 
separable state of n qubits. The choice of U$ is optimized 
in order to maximize the success probability of the search 
process for the given state p. 



functions of the form of G(ip) satisfy the conditions 
above for an entanglement monotone. More specifically, 
there exist two specific purifications, |<^o) and \i()q), f° r 
which F(p,a) = | <<^o |-0o) | 2 - Thus, G(p) = if and 
only if F(p,a) = 1. In this case K^olV'o)! 2 = 1, or 
\4>o) = e^l^o); thus p = a. Since a is separable then 
so is p, and the first condition is satisfied. 

In order to prove the second condition one can use the 
monotonous quality of the fidelity under trace preserving 
operations Q . This means that for every complete set of 
operators 



in) 



(19) 



where and {|?q)} are orthonormal basis states and 

Ur and Uq are unitary operators, in R and Q, respec- 
tively. 

For any separable mixed state, a of n qubits, the uni- 
tary operators Ur and Uq and the operator y/a provide 
a convenient parametrization of the pure state \<j>). This 
follows from the fact that Eq. (JIJ can also be expressed 



as a = 



where a,,, 



k = 1, . 



is the 



density operator of a mixed state of the kth qubit. To 
obtain the operator y/a, one needs to diagonalize a, ac- 
cording to VaV^ — D, where V is a suitable unitary 
operator, resulting in a diagonal matrix D. Taking its 
square root one obtains another diagonal matrix d with 
matrix elements da = y 'Da. The process is completed 
by \[& = V^dV. The operator s/a is not a required to 
be unitary and is not used in the quantum circuit. 

In order to complete the construction of the quantum 
circuit, one needs the operator ul that transforms \rj) 
into a purification \<p) of a separable state a. To this end 
we construct two bases of the space spanned by 2n qubits. 
The first basis is obtained from the computational basis, 
\i) by appling the Hadamard transform on all qubits, 
namely \rji) = H® 2n \i). Note that 1 770) = |»7). The second 
basis, \4>i), can be constructed using the Gram-Schmidt 
algorithm, starting with 1 0o) = 10)- The unitary operator 
U\ = J^i transforms the state \rj) into the state 

\4>), which is a purification of a separable state, cr, of n 
qubits. 



F(p, a) < F MipM},^ MiaMjj , (17) 

where the separable state a remains separable under the 
transformation. As a result, 



maxF(p, cr) < mfrxf f ^ M.pM^aj . (18) 

Therefore, G(p) cannot increase under such transforma- 
tions. 



B. The Unitary Operator U$ and its Operational 
Interpretation 

Consider a mixed state a of n qubits. The state a 
can be purified into a pure state of 2n qubits, half of 
them associated with the original subspace Q and the 
rest with the added subspace R |^y. Such purification 



V. SUMMARY AND DISCUSSION 

The Groverian entanglement measure for multiple 
qubits, previously introduced for pure states, was gen- 
eralized to the case of mixed quantum states. This gen- 
eralization provides an operational measure of entangle- 
ment for both pure and mixed states. The operational 
interpretation is based on the fact that the Groverian 
measure of a state p is related to the success probability 
of Grover's search algorithm when p is used as the initial 
state, following a suitable pre-processing stage. When p 
is a mixed state, the pre-processing includes a purifica- 
tion procedure followed by the application of a certain 
unitary operator. In case that the given state p is a pure 
state, G(p) coincides with the Groverian measure for pure 
states, introduced in Ref. [Tsj . 

It might seem surprising that in order to evaluate the 
entanglement of an rt-qubit system one needs a search 
space of In qubits. This can be explained by considering 
mixed states as open systems, where, in some cases, the 
mixture represents an effective entanglement with exter- 
nal qubits. Purification of the mixed state p to 2n qubits 
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enables us to create a closed system with no entangle- 
ment to any external qubits, in which all the relevant 
information is maintained. The Groverian measure for 
mixed states enables to fully capture the entanglement in 
p. Mathematically, this measure satisfies the conditions 
required from an entanglement measure. In particular, 
it vanishes for any separable mixed state of n qubits, it 
is invariant under local unitaries and monotonically de- 
creasing (in the weak sense) under LOCC. 

Quantum algorithms are designed to start with a well 
defined initial state. The final state, just before the mea- 
surement is taken, can be either a basis state or a super- 
position of basis states. For example, in Grover's algo- 
rithm with a single marked state, the desired final state 
(namely the marked state) is a basis state. In Grover's 
algorithm with several marked states, as well as in Shor's 
factoring algorithm, the desired final state is a superpo- 
sition. The analysis presented in this paper can be gen- 
eralized by replacing Grover's algorithm by some other 
quantum algorithm. If in the replacement algorithm the 
desired final state is a basis state, G(p) will not depend on 
the specific algorithm. However, for algorithms such as 
Shor's algorithm, the maximal success probability may 
not coincide with the expression used in the Groverian 
measure. Yet, in the special case of Grover's algorithm 
with several marked states, the Groverian measure still 
holds H3. 

Recently, the Groverian measure was applied in order 
to evaluate the entanglement in certain pure quantum 
states of multiple qubits 20] . A convenient parametriza- 
tion was developed that enables analytical calculations 
of G(ip) for some pure states of high symmetry. In order 
to evaluate it for arbitrary states, a numerical minimiza- 
tion procedure, based on the steepest descent algorithm 
was developed. Using this procedure, the entanglement 
of intermediate states, generated during the evolution of 
Grover's algorithm, was calculated. It was found that 
even if the initial state and the target state are product 
states, in intermediate stages of the algorithm, highly 
entangled states are generated, in agreement with ear- 



lier studies in which other measures were used [28|, |29j . 
This result is interesting in the context of attempts to 
examine the role of entanglement in quantum algorithms 
and specifically in Grover's search algorithm. In partic- 
ular, recent studies have shown that an implementation 
of Grover's algorithm using classical media, namely, in 
which quantum entanglement docs not play a role, would 
require an exponentially larger overhead compared to the 
quantum case [13, El ■ 

Unlike the case of pure states of two qubits, in which 
the von Neumann entropy provides a complete character- 
ization of the entanglement, multiple qubit states support 
a large number of different measures [2^, |23, yj, |33|, y4| • 
It seems that the issue of what measure is relevant de- 
pends on the specific physical or operational context in 
which it is used. In particular, the Groverian entangle- 
ment measure is motivated by a quantum algorithm. It 
thus appears to be a suitable measure for the evaluation 
of the entanglement that is produced during the evolu- 
tion of quantum algorithms. The actual evaluation of 
entanglement measures turns out to be a difficult com- 
putational problem. This is due to the fact that these 
measures are typically defined as an extremum of some 
multi-variable function. A singular result in this context 
is the explicit formula for the entanglement of formation 
of mixed states of two qubits, obtained in Refs. [T^.ll3| . 

A related operational interpretation of the fidelity, 
which is also based on Uhlmann's theorem, was intro- 
duced in Ref. [3{|. In that case the fidelity F(p, a) of the 
output states of a noisy channel provides an upper bound 
on the overlap of the input states, under the assumption 
that they were pure. 

The generalization of the Groverian measure to mixed 
states may provide further insight into the role of en- 
tanglement in making quantum algorithms powerful. It 
would be interesting to use this measure to evaluate the 
entanglement generated by quantum algorithms using 
mixed states, particularly when decoherence effects are 
taken into account. 
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